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Abstract 



It is shown that the i?-parity violating decays of the lighter top squarks (ti) trig- 



gered by the lepton number violating couplings A^ 33 , where the lepton family index 
\ i = 1-3, can be observed at the LHC via the dilepton di-jet channel even if the 

coupling is as small as 10~ 4 or 10~ 5 , which is the case in several models of neu- 
trino mass, provided it is the next lightest supersymmetric particle(NLSP) the lightest 
neutralino being the lightest supersymmetric particle(LSP). We have first obtained a 
fairly model independent estimate of the minimum observable value of the parameter 
(Pij = BR(t—* ifb) x BR(t* -> Ijb)) at the LHC for an integrated luminosity of lOfb" 1 
as a function of m~ by a standard Pythia based analysis. We have then computed 
the parameter P« in several representative models of neutrino mass constrained by the 
neutrino oscillation data and have found that the theoretical predictions are above the 
estimated minimum observable levels for a wide region of the parameter space. 
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1 Introduction 



Neutrino oscillation experiments [1] have confirmed that neutrinos indeed have very tiny 
masses, several orders of magnitude smaller than any other fermion mass in the Standard 
Model (SM). The tiny masses of the neutrinos, however small, provide evidences of new 
physics beyond the SM. 

Neutrinos can be either Dirac fermions or Majorana fermions depending upon whether 
the theory is lepton number conserving or violating. In the SM, as originally proposed by 
Glashow, Salam and Weinberg, neutrinos are massless since right handed neutrinos and 
lepton number violating terms are not included. 

Both i?-parity conserving (RPC) or .R-parity violating (RPV) minimal supersymmetric 
extension of the SM (MSSM)[2] are attractive examples of physics beyond the SM. In general 
the MSSM may contain RPC as well as RPV couplings. The latter include both lepton 
number and baryon number violating terms which result in catastrophic proton decays. One 
escape route is to impose i?-parity as a symmetry which eliminates all RPV couplings. This 
model is generally referred to as the RPC MSSM. However, neutrino masses can be naturally 
introduced in this model only if it is embeded in a grand unified theory (GUT) [3]. Tiny 
Majorana neutrino masses are then generated by the see-saw mechanism [4]. Proton decay 
is a crucial test for most of the models belonging to this type. 

However, an attractive alternative for generating Majorana masses of the neutrinos with- 
out allowing proton decay is to impose a discrete symmetry which eliminates baryon number 
violating couplings from the RPV sector of the MSSM but retains the lepton number violat- 
ing ones. The observation of neutrinoless double beta decay[5] and absence of proton decay 
may be the hallmark of such RPV models of neutrino mass. 

The GUT based models though very elegant have hardly any unambiguous prediction 
which may tested at the large hadron collider (LHC). In contrast the RPV models of neutrino 
mass are based on TeV scale physics and, consequently, have many novel collider signatures. 

The observables in the neutrino sector not only depend on the RPV parameters but also 
on the RPC ones like the masses of the superpartners generically called sparticles. Thus 
the precise determination of the neutrino masses and mixing angles in neutrino oscillation 
experiments together with the measurement of sparticle masses and branching ratios (BRs) at 
collider experiments can indeed test the viability of the RPV models quantitatively. Moreover 
the collider signatures of this model are quite distinct from that of the RPC model. In this 
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paper our focus will be on a novel signature of a RPV model of v mass which can be easily 
probed at the early stages of the upcoming LHC experiments. 

In the RPC models the lightest supersymmetric particle (LSP) decays into lepton number 
violating channels producing signals with high multiplicity but without much missing energy 
which are in sharp contrast with the signals in a typical RPC model. In RPV MSSM the 
sparticles other than the LSP can also directly decay via lepton number violating channels 
which may lead to spectacular collider signatures. However, in a typical model of neutrino 
mass consistent with the oscillation data such couplings turn out to be so small [6] that the 
RPC decay of the sparticles overwhelm the RPV decays. Thus the LSP decay is the only 
signature of i?-parity violation. 

However, the scenario changes dramatically if we consider the direct RPV decay of 
the lighter top squark (t\) [7, 8, 9, 10] with the assumption that ti it is the next lightest 
supersymmetric particle(NLSP) while the lightest neutralino (xT) i s the LSP. The theoretical 
motivation for the ti-NLSP scenario is the fact that it's superpartner - the top quark- is 
much heavier than any other matter particle in SM. This large top mass (m*) leads to a 
spectacular mixing effect in the top squark mass matrix which suppresses the mass of the 
lighter eigenstate [2]. We assume that tx-NLSP decays via the loop induced mode ti — > c\\ 
[11] and the four body [12] decay mode, which occurs only in higher order of perturbation 
theory, with significant BR. The validity of this assumption will be justified later. The RPV 
decays can now naturally compete with the RPC ones in spite of the fact that couplings 
underlying the former modes are highly suppressed by the v oscillation data [13]. 

The lighter top squark decays into a lepton and a 6-jet via RPV couplings A' 33 are listed 
below: 

a) t x l+b ; b) t\ l~b (1) 

where 2=1-3 corresponds to e, /i and r respectively. Our signal consists of opposite sign 
dileptons(OSDL), two hard jets with very little pr- These modes dominate, e.g., in many 
RPV models where neutrino masses are generated at the one loop level by the A^ 33 cou- 
plings, where i is the lepton index and 3 stands for quarks or squarks belonging to the third 
generation (see below). 

We take the lowest order QCD cross section of top squark pair production which depends 
on mqr only. Requiring that the significance of the signal over the SM background be at least 
5 a level for an integrated luminosity of 10 fb~ l or smaller, we can then put fairly model 
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independent lower limits on the products of the BRs (PBRs) of the RPV decay modes in 
Eq. 1. In our analysis both the signal and the backgrounds are simulated with Pythia. As 
expected the range of vnq which can be probed at the LHC is significantly larger compared 
to the reach of Tevatron RUN II[14, 15]. The details of our simulations will be presented in 
the next section. 

In principle the viability of probing any RPV model of neutrino mass with the above 
characteristics at the LHC can be checked by computing PBRs in respective models, and 
comparing with the estimated lower limits. For the purpose of illustration we have considered 
in section 3 a model based on three bilinear RPV couplings (/ij) and three trilinear couplings 
(A^ 33 ) at the weak scale [16] and have carried out the above check. It is gratifying to note 
that most of the parameter space allowed by the neutrino oscillation data can be probed 
by the early LHC experiments with an integrated luminosity of 10 fb _1 (see section 3). 
Moreover, the constraints from oscillation data indicate that the A^ 33 couplings should have 
certain hierarchical pattern leading to distinct collider signatures [17] . This hierarchy among 
the couplings can be qualitatively tested by observing the relative sizes of signals involving 
different OSDL signals. 

The summary, the conclusions and future outlooks are in the last section. 

2 The signals and the SM backgrounds 

The production and decay of the lighter top squark pairs are simulated by Pythia[18]. Initial 
and final state radiation, decay, hadronization, fragmentation and jet formation are imple- 
mented following the standard procedures in Pythia. We have considered only the RPV 
decay modes of t\ via the couplings A^ 33 ,i = 1-3 (Eq. 1) and in this section their BRs are 
taken to be free parameters. We have used the toy calorimeter simulation (PYCELL) in 
Pythia with the following criteria: 

• The calorimeter coverage is \rj\ < 4.5. The segmentation is given by At] x A<p = 
0.09 x 0.09 which resembles a generic LHC detector. 

• A cone algorithm with AR(j,j) = \J Aif + A0 2 = 0.5 has been used for jet finding. 

• Jets are ordered in Et and E^ min = 30GeV. 
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Various combinations of OSLDs in the final state are selected as follows: 

• Only tau leptons decaying into hadrons are selected provided the resulting jet has 
P T > 30 GeV and \r]\ < 3.0. 

• Leptons (/ = e, /i) are selected with Pt > 20 GeV and \rj\ < 2.5. 

The following selection criteria(SC) are used for background rejection : 

• The r-jets are tagged according to the tagging efficiencies provided by the CMS 
collaboration [19] (Fig. 12.9) (SCI). Hadronic BR of the r is also included in the cor- 
responding efficiency. For e and fi SCI is the lepton-jet isolation cut. We require 
AR(l,j) > 0.5. The detection efficiency of the leptons are assumed to be approximately 
100% for simplicity. 

• Events with two isolated leptonic objects (e,/i or tagged r-jets ) are rejected if Pt < 150 
GeV, where I = e or fi (SC2) or (r) <100 GeV, where £^ (r) is the E T of the r jet. 

• We select events with exactly two jets other than the tagged r-jets (SC3). The event 
is rejected if the additional jets have Pt < 100 GeV (SC 4). 3 

• Events with missing transverse energy (Et) > 60 GeV are rejected (SC5). 

Through SCI we have severely constrained the transverse momentum of two leptons / = 
e, \i to reject the leptons coming from the leptonic decays of the tau. Moreover such a strong 
cut reduces most of the SM backgrounds significantly. We have considered backgrounds 
from: WW, WZ, ZZ, ti, Drell-Yan (DY) and QCD events. The missing energy veto plays 
a crucial role to tame down WW and ti backgrounds as they are rich in missing energy. 
Mistagging of light jets as r-jets is a major source of background to di-tau events. We have 
taken this into account. However, if we also employ fe-tagging then this background can be 
brought under control to some extent. 

In our work b tagging has been implemented according to the following prescription. A 
jet with \r)\ < 2.5 matching with a B-hadron of decay length > 0.9 mm has been marked 
tagged. The above criteria ensures that e& ~ 0.5 in ti events, where e& is the single 6-jet 

3 It is expected that this cut would also suppress the SUSY backgrounds due to, e.g., q, g production. 
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tagging efficiency (i.e., the ratio of the number of tagged 6-jets and the number of taggable 
6-jets in tt events). 

The leading order(LO) cross-sections for t\ — 1\ pairs presented in Table 1 are computed 
using calcHEP (version 2.3.7) [20]. 



Signal 


240 


300 


400 


450 


500 


a (jab) 


14.6 


4.8 


1.1 


0.58 


0.32 



Table 1: t\ - t* pair production cross section (a) at the LHC for different . 

In Table 2 we have presented the combined efficiencies of SCI - SC5 in steps. The first 
column of Table 2 shows signals with different topology of final states. Here e e X, ttX, erX 
and e/j,X represent final states without b-jet tagging. The cumulative efficiency of each SC 
for = 400 GeV is presented in the next five columns. However, we have not separately 
presented the efficiencies corresponding to final states with muons as we have assumed that 
both e and \i are detected with approximately 100% efficiency. 

Table 3 contains the effect of b-jet tagging on different final states. We have used 
the notations 0b, lb and 26 to specify signal events with zero, one and two tagged b-jets 
respectively. From this Table it is also evident that the efficiencies increase for larger vnq 
since the Pt cut on leptons become less severe. This compensates the fall of the cross section 
with increasing vrtq to some extent. 



t x t\ 


ei 


£2 


£3 


£4 


C5 


eeX 


0.93708 


0.292239 


0.087228 


0.043344 


0.032823 


ttX 


0.251343 


0.111546 


0.033201 


0.031554 


0.008955 


e\iX 


0.94101 


0.295239 


0.088216 


0.043415 


0.033060 


erX 


0.474948 


0.180945 


0.053820 


0.044793 


0.016965 



Table 2: Efficiency table for m Ti = 400 GeV. 

In Table 4 we have shown the effect of cuts on the background from tt events. SC2 is very 
effective in reducing this background significantly. Moreover this background is accompanied 
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SIGNAL 


mr ( GeV) 
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300 


400 


450 


500 


eeOb 


0.00032 


0.00066 


0.00189 


0.00234 


0.00255 


e e lb 


0.00121 


0.00330 


0.01116 


0.01461 


0.01580 


e e 2b 


0.00176 


0.00509 


0.01984 


0.02620 


0.03121 


eeX 


0.00328 


0.00905 


0.03282 


0.04315 


0.04957 


T T Ob 


0.00059 


0.00073 


0.00112 


0.00091 


0.00097 


r r lb 


0.00153 


0.00284 


0.00363 


0.00351 


0.00391 


r t 2b 


0.00126 


0.00226 


0.00421 


0.00450 


0.00522 


TT X 


0.00338 


0.00582 


0.00896 


0.00892 


0.01098 


t e Ob 


0.00045 


0.00081 


0.00148 


0.00142 


0.00142 


t e lb 


0.00135 


0.00307 


0.00667 


0.00705 


0.00717 


t e 2b 


0.00126 


0.00346 


0.00882 


0.00997 


0.01078 


t e X 


0.00308 


0.00734 


0.01697 


0.01843 


0.01936 


H e Ob 


0.000315 


0.00067 


0.00190 


0.00235 


0.00257 


H e lb 


0.00123 


0.00334 


0.01125 


0.01469 


0.01635 


H e 2b 


0.00178 


0.00512 


0.01992 


0.02625 


0.03129 


fjb e X 


0.00332 


0.00912 


0.03306 


0.04329 


0.05021 



Table 3: Final efficiencies for different (including b-tagging if implemented). 
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1.88 x 10" 3 


9.3 x 10" 6 


2.95 x 10" 5 


9.50 x 10~ 7 


2.00 x 10" 7 



Table 4: Efficiency table for tt process 

by large amount of and SC5 also reduces it significantly. Since tt decays contain two b 
quarks, b- tagging is not very effective here and has not been included in Table 4. 



6 



QCD 


ei 


£2 


£3 


64 


£5 


ee 


1.16 x 10" 5 














TT 


9.10 x 10" 3 


4.02 x 10" 3 


1.05 x 10" 3 


2.85 x 10" 4 


2.10 x 10" 4 


e/i 
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Table 5: Efficiency table for the QCD process in the Pt bin: 400 GeV < Pt < 1000 GeV. 

Table 5 presents another important background arising from the 2^2 processes due to 
pure QCD interactions for 400 GeV < pT < 1000 GeV, where Pt is the transverse momentum 
of the two partons in the final state . However, SC2 completely kills all backgounds except 
for those with the di-r final states. The latter background, mainly due to mistagging of light 
flavour jets as r-jets, affect the di-r signal very seriously . The mistagging probability has 
also been taken from [19] (Fig. 12.9). 

This background is very large, as expected, since the QCD cross-section is very large. 
The leading order cross-sections have been computed by Pythia in two px bins : (i) 400 GeV 
< p T < 1000 GeV and (ii) 1000 GeV< p T < 2000 GeV . We have chosen the QCD scale 
to be v£- The corresponding cross-sections being 2090 pb and 10 pb respectively. Beyond 
2000 GeV the number of events are negligible. We shall discuss later how the visibility of 
the di-r signal can be improved by employing b tagging. 

In Table 6 we have computed the numerically significant backgrounds of all types for C = 
10 fb _1 . Here '-' denotes a vanishingly small background. It is clear from this table that only 
ti and QCD backgrounds are relevant. The LO cross-sections in the second row of Table 6 
except for the QCD processes have been computed using calcHEP (version 2.3.7) [20]. Due to 
very strong cut on Pt of highest two leptons SC2 DY type backgrounds become vanishingly 
small. Moreover, SC3 and SC4 finally reduce it to zero. Other backgrounds like WW, WZ 
and ZZ become vanishingly small mainly due to SC2. 

The Product Branching Ratio (PBR) is defined as : 

P tj = BRfo -> Ifb) x BR(t\ -> IJb) (2) 

where % or j can run from 1-3 corresponding to e, fi and r respectively. The Minimum 
Observable Product Branching Ratio(MOPBR = P™ n ) corresponds to S/ y(B) > 5, where 
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0.80 
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- 


r r 26 
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0.20 
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0.40 






fi e lb 








0.80 






fi e 2b 
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fi e 
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2.00 







Table 6: Total number of all types of backgrounds survived after all cuts . 

S and B are the number of signal and background events respectively. However, for a typical 
signal with negligible background we have required S > 10 as the limit of observability and 
MOPBR is computed accordingly. 

For a given L the MOPBR for each process is computed from Table 3 and Table 6 by 
following expression: 
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where P^ is already defined in Eq. 2. a b and e b (not to be confused with e&, the 6-jet tagging 
efficiency) denote the cross section and the efficiency of background of type b . Similarly e 
is the final efficiency for the signal, rj is 2 for i ^ j and t] is 1 for i — j. The integrated 
luminosity £ is taken to be 10 fb _1 . The estimated MOPBRs are given in Table 7(without 
6-jet tagging) and Table 8 ( with two tagged 6-jets). 

We remind the reader that in Table 7 and Table 8 a signal is assumed to be observable 
if S > 10 even if B is < 4. In Table 7 and Table 8 a 'x' indicates that corresponding channel 
can not be probed. 

Our conclusions so far have been based on LO cross sections. If the next to leading order 
corrections are included the t\ — t\ production cross section is enhanced by 30 - 40 % due 
to a K-factor[21]. It is then clear from Eq. 3, that the estimated MOPBR would remain 
unaltered even if all significant background cross sections are enhanced by a factor of two 
due to higher order corrections. 
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X 


X 
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1.0 


1.1 


1.4 


2.0 


3.4 


p^ in {%) 


1.1 


1.4 


2.7 


4.7 


8.0 



Table 7: Minimum value of PBR estimated from the sample without b tagging . 

We present in Fig. 1 the distribution (unnormalised) of invariant mass of a electron -jet 
pair in the dielectron-dijet sample without b-tagging for mj = 300 GeV. We first reconstruct 
invariant mass for all possible electron-jet pair. Among these pairs We select the two such 
that the difference in their invariant mass is minimum. We then plot the higher of the two 
invariant masses. This peak, if observed, would unambiguously establish the lepton number 
violating nature of the underlying decay. In contrast if neutralino decay is the only signal of 
P-parity violation, then this information may not be available. For example, if Xi ~ * u bb is 
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m Ti (GeV) 
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4.1 
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6.6 


10.0 


Pg»{%) 


9.4 


16.0 


37.5 


66.4 


X 


P% in (%) 


1.9 


2.0 


2.3 


3.3 


5.0 




2.7 


3.0 


5.2 


8.6 


14.5 



Table 8: Minimum value of PBR estimated from the 2-6 tagged sample. 
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Figure 1: The invariant mass distribution for m^r = 300 GeV. 

the dominant decay mode of the LSP via the A -33 coupling then the lepton number violat 
nature of the decay dynamics will be hard to establish. 
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In the next section we shall calculate the PBR for different signals in a realistic models 
of neutrino mass constrained by the neutrino oscillation data and examine whether the 
predictions exceed the corresponding MOPBR estimated in this section. Our main aim is to 
illustrate that the LHC experiments will be sufficiently sensitive to probe these models and 
not to make an exhaustive study of all possible models. 

3 Model Calculations 

The collider signatures considered in the last section arise only in models with non-vanishing 
trilinear A^ 33 type couplings at the weak scale. However, consistency with neutrino oscilla- 
tion data require the introduction of more RPV parameters (bilinear superpotential terms, 
bilinear soft breaking terms etc) [22]. In fact the list of possible choices is quite long. It is 
expected that the constraints on the A' couplings in the most general model imposed by the v 
- oscillation data will be considerably weaker and the observability of the resulting dilepton- 
dijet signal will improve. Thus we have restricted ourselves to models with a minimal set 
of parameters capable of explaining the oscillation data with rather stringent constraints on 
the A' couplings. 

We work in a basis where the sneutrino vevs are zero. It is assumed that in this basis only 
three nonzero bilinear (/Xj) and three trilinear (A^ 33 ) couplings, all defined at the weak scale, are 
numerically significant. In this framework the neutrino mass matrix receives contributions 
both at the tree and one loop level. It should be emphasised that the tree level mass matrix, 
which is independent of Aj33 couplings, yields only two massless neutrinos. Thus the interplay 
of the tree level and one loop mass matrices is essential for consistency with the oscillation 
data. 

The char gino- charge lepton, the neutralino - neutrino and other relevant mixing matrices 
in this basis may be found in [22]. In principle the diagonalization of these matrices may 
induce additional lepton number violating couplings which can affect the BRs of the top 
squark decays considered in this paper. For example, the RPC coupling t\ — t — W3 may 
induce new RPV vertices through W3 — v mixing. However, it was shown in [23] that the 
new modes induced in this way would have negligible BRs. As a result the approximation 
that the decays of the top squark NLSP are driven by the A^ 33 couplings only is justified. 

In addition to the RPV parameters the neutrino masses and mixing angles depends on 
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RPC parameters. In this paper we shall use the following popular assumptions to reduce 
the number of free parameters in the RPC sector: i) At the weak scale the soft breaking 
mass squared parameters of the L and R-type squarks belonging to the third generation are 
assumed to be the same( the other squark masses are not relevant for computing neutrino 
masses and mixing angles in this model), ii) We shall also use the relation M 2 ~ 2 Mi at 
the weak scale as is the case in models with a unified gaugino mass at Mq. Here M\ and 
M 2 are respectively the soft breaking masses of the U(l) and SU(2) gauginos respectively. 

The tree level neutrino mass matrix and, hence, the predicted neutrino masses depends 
on the parameters of the gaugino sector (through the parameter C[16, 17]). They are M 2 , 
Mi , fi (the higgsino mass parameter) and tan (3 = v 2 /vi, where V\ and v 2 are the vacuum 
expectation values (vevs) for the down type and the up type neutral higgs bosons respectively. 
We remind the reader that for relatively large tan j3s the loop decay overwhelms the RPV 
decay [10, 12]. We have, therefore, restricted ourselves to tan/3 = 5-8. It is also convenient 
to classify various models of the RPC sector according to the relative magnitude of M 2 and 
/1. If Mi < M2 <C /i, then the lighter chargino (xt )> the LSP (xi) an d the second lightest 
neutralino (x 2 ) are dominantly gauginos. Such models are referred to as the gaugino-like 
model. On the other hand in the mixed model (Mi< M 2 ~ //), xt an d X 2 are admixtures 
of gauginos and higgsinos. In both the cases, however, x? is purely a bino to a very good 
approximation. There are models with Mi,M2 ^> /i in which Xi j Xi an d X2 are higgsino 
- like and all have approximately the same mass (~ fi). It is difficult to accommodate the 
top squark NLSP in such models without fine adjustments of the parameters. Thus the 
LSP decay seems to be the only viable collider signature. One can also construct models 
wino or higgsino dominated LSPs. However, the ti-NLSP scenario cannot be naturally 
accommodated in these frameworks for reasons similar to the one in the last paragraph. 

The one loop mass matrix, on the other hand, depends on the sbottom sector (through 
the parameter K 2 [16, 17]). This parameter decreases for higher values of the common squark 
mass for the third generation. From the structure of the mass matrix it then appears that 
for fixed C, identical neutrino masses and mixing angles can be obtained for higher values 
of the trilinear couplings if K 2 is decreased. Thus at the first sight it seems that arbitrarily 
large width of the RPV decays may be accommodated for any given neutrino data. This, 
however, is not correct because of the complicated dependence of the RPV and loop decay 
BRs of ti on the RPC parameters and certain theoretical constraints. The common squark 
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mass cannot be increased arbitrarily without violating the top squark NLSP condition. Of 
course larger values of the trilinear soft breaking term A t may restore the NLSP condition. 
But larger values of A t tend to develop a charge colour breaking( CCB ) minimum of the 
scalar potential [24]. Finally the pseudo scalar higgs mass parameter Ma can be increased 
to satisfy the CCB condition. But as noted earlier [17] that would enhance the loop decay 
width as well and suppress the BRs of the RPV decay modes. 

We have chosen the following RPC scenarios : A) The gaugino dominated model and 
B)The mixed type model. The choice of RPC parameters for model A) and model B) are: 
A) Mi = 195.0, M 2 = 370.0, // = 710.0, tan/5 = 6.0, A t = 1100.0, A b = 
1000.0, Mq(common squark mass )=450.0, Mj (common slepton mass ) = 400.0 and 
M A = 500.0 and B) M x = 170.0, M 2 = 330.0, n = 320.0, tan/3 = 6.0, A t = 1045.0, A b = 
1000.0, Mg=450.0, Mj = 400.0 and Ma = 200.0, where all masses and mass parameters are 
in GeV . Both the scenarios correspond to = 240 GeV and ti is the NLSP. It should be 
noted that the slepton mass is specified to ensure that the ti is the NLSP. It does not affect 
the neutrino mass matrix. 

Even if ti is the NLSP the following modes may compete with the RPV decays and 
overwhelm it: 

a) h^txl; b) h^bWxl c)ti-cx?; d) h -> ffbft (4) 

In the parameter spaces we have worked with the mode a) is kinematically disallowed. 
The second mode is highly suppressed if the LSP is Bino dominated as is assumed in this 
analysis. Thus in the scenario under consideration only modes c) and d) may compete with 
RPV decays of t\. In this section we have computed the PBRs taking into account the 
competition among the above three modes. 

Next we have randomly generated bilinear and trilinear RPV couplings, /ij and A^ 33 . 
Then these parameters are constrained by v oscillation data which allows very few sets of 
RPV parameters for the above RPC parameters. Most of the allowed trilinear RPV couplings 
lie within 10~ 4 — 10~ 5 . Finally the relevant PBRs have been calculated in model (A) and 
(B). 

In Table 9 we present several representative sets of trilinear RPV parameters allowed 
by v oscillation data and the corresponding PBR. A '-' indicates that the predicted PBR is 
negligible. As noted before two of the couplings turns out to be large while the third one 
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is suppressed due to oscillation constraints. It turns out that the PBR's involving the large 
couplings are larger than the corresponding MOPBRs estimated in the last section even 
without b-tagging(see Table 7). The only exception is P33 which cannot be probed without 
b-jet tagging (see Table 8) 



A'i 33 [xl0- 5 ] 


A' 23 3[xl0- 5 ] 


A' 3 33[X10- 5 ] 


P11 


P22 


-P33 


P12 


P23 


Pl3 


Model A 




1.6 


8.3 


10.0 




5.1 


10.8 


0.2 


7.4 


0.3 


7.5 


0.7 


9.2 


4.9 




11.4 




0.1 


7.5 


4.6 


4.5 


0.3 


6.8 


6.8 




6.7 






Model B 




11.9 


0.99 


15.0 


4.2 




10.6 






6.6 


0.59 


13.6 


16.8 




4.3 


10.2 




6.6 




7.3 


7.4 


0.9 


6.3 


6.6 




6.4 


0.1 


0.1 



Table 9: Trilinear RPV couplings allowed by v oscillation data and the corresponding PBRs 
computed in models A and B (see text) with mj = 240 GeV . 

For larger , there exists allowed RPV parameter space with observable PBRs at the 
early LHC runs. However, if we go beyond mqr = 500 GeV the di-tau channel cannot be 
probed even with b -tagging. Nevertheless, observation of the e — r and the \x — r channel 
will provide evidence for a relatively large A333. We present in Table 10 for mq = 500 GeV. 
The RPC parameters corresponding to a Gaugino model are chosen to be: 
Mi = 475.0, M 2 = 860.0, li = 1650.0, tan/3 = 6.0, A t = 995.0, A b = 1000.0, M g -=575.0, 
Mi = 525.0 and Ma = 300.0, where all masses and mass parameters are in GeV. 



A'i3 3 [xl0- 5 ] 


A' 233 [xl0- 5 ] 


A' 3 33[xl0- 5 ] 


Pn 


P22 


-P33 


P12 


P23 


Pi* 


9.1 


4.0 


6.4 


20.7 




5.1 


4.0 


2.0 


10.3 


4.4 


10.9 


5.6 




31.4 


2.2 


5.1 


8.3 


1.3 



Table 10: Same as Table 9 for = 500 GeV. 
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We have checked that even for mq > 500 GeV there exits RPV parameter space allowed by 
oscillation data which leads to observable dilepton -dijet signals in early LHC experiments. 



/ \ / \max 
l A 133J 


1 \ I \max 
\ A 233J 


( \ l \ max 
l A 333J 


92 


2176 


2086 


Pn, P13 


56 


-P22, -P23 


1376 


P33, -P23 


874 


Pn, P12 


15 


-P22, P12 


119 


-^33, P\3 


27 


Pn 


2 


P22 


664 


P22, P23 


274 


P12 


10 


** 


17 


Pn, P13 


25 


** 


9 






P23 


45 










P33 


304 










** 


537 



Table 11: Number of allowed solutions in the mixed model (m-^ = 240GeV) consistent with v 
oscillation data which satisfy the MOPBR given in Table 7 and Table 8. The above numbers 
are estimated for C = 10fb _1 . 

We have randomly generated 10 9 sets of RPV parameters in the mixed model with 
mj = 2A0GeV. Out of these only 4354 are consistent with the ^-oscillation data. These 
solutions can be further classified into three groups according to the highest value of A^ 33 . 
The three columns in Table 11 correspond to these groups. The first column in Table 11 
contains detailed information about the flavour structure of the RPV couplings in the 92 
solutions with the hierarchy A' 133 > A' 233 ,A 333 . The next few rows display the number of 
solutions with predicted PBRs in different channels above the observable limits as given in 
Table 8. For example, the third row indicates that signals in ee + 2jets and er + 2jets 
channels are observable with 10 fb _1 of data in 56 solutions. These channels, if observed, 
would further reveal that A' 133 > A 333 > A' 233 . On the other hand observable signals in 
ee + 2jets and e/x + 2jets channels as given in the next row would indicate the hierarchy 

A 133 ^ A 233 ^ A 233' 

If only one channel, say the ee + 2 jets, is observed one can conclude that A' 133 >> 
A' 233 , A 333 (see row 5). On the other the observation of the e/x + 2 jets signal only (see row 6) 
would indicate A' 133 ~ A' 233 >> A 233 . The channel e/i + 2 jets dominates over the ee + 2 jets 
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or the + 2jets channel because of the factor of two which enhances the number of events 
when leptons of two different flavours with all possible charge combinations are observed. 
Finally the seventh row with '**' indicates that no signal can be observed with £ = 10 fb _1 
The information in the next two columns are presented following the format and similar 
inferences about the hierarchy of the A^ 33 can be drawn from the lepton flavour content of 
the final states. We have verified that for £ = lOOftT 1 all solutions well predict atleast one 
Pij above the corresponding P™ n . 

4 Conclusion 

In conclusion we reiterate that the OSDL signals with same or different flavours of leptons 
(e, \i or tau-jets) plus two additional jets arising from RPV decays of t\ - t\ pairs produced 
at the LHC would be a promising channel for probing the RPV couplingA- 33 (see Eq. 1 and 
the discussions following it). This is true in general if t\ happens to be the NLSP, which is 
a theoretically well motivated scenario. This signal is especially interesting in the context of 
RPV models of neutrino mass. A part of our analysis (section 2), however, is fairly model 
independent since the size of the signal is necessarily controlled by the production cross 
section of the t\ - t* pair as given by QCD and the product branching ratio (see Eq. 2). 
The model independent estimates of P^ m (see Eq. 3) corresponding to observable signals 
for different m^s (see Eq. 3) for an integrated luminosity of 10 fb^ 1 are presented in Table 
7 and Table 8 using the Monte Carlo event generator Pythia. We have optimized the cuts 
for m-jr = 240 GeV. However, for even larger values of vrtq the signal efficiencies increase for 
the same set of cuts keeping the background events almost negligible. Top squark masses 
in the vicinity of 500 GeV yield observable signals in this channel for realistic models of 
m u . Although our calculations are based on LO top squark pair production cross sections 
we emphasise that the inclusion of NLO corrections are likely to yield even larger estimates 
of P™ n as argued in section 2. 

We have further noted that inspite of the combinatorial backgrounds, the invariant mass 
distribution of the lepton (e or /i)-jet pair shows a peak at mj (see Fig.l). This peak, 
if discovered, will clearly establish the lepton number violating nature of the underlying 
interaction. This may not be possible if neutralino decays happen to be the only RPV 
signal. 
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In models of z/-mass, the underlying A' couplings turn out to be very small. If A^ 33 
contributes to the one loop z/-mass matrix, it is typically of the order of 10~ 4 - 10~ 5 due 
to constraints imposed by the //-oscillation data. Even if A' is so small the RPV decay of 
the ti-NLSP may have sizable BRs over a large region of the parameter space because the 
competing loop induced decay (Eq. 4c) or the four body decay (Eq. 4d) of t\ also have 
suppressed widths. For the purpose of illustration we have considered a specific model of 
z/-mass[16] with parameters constrained by the //-oscillation data. It is interesting to note 
that in this model most of the theoretically predicted -Pj/s (Eq. 2) for several representative 
choices of RPC parameters turn out to be larger than the P^ m 's estimated in section 2 
for £= 10 fb~ l . For larger C almost all solutions yield Pi/s at the observable level. The 
relative size of the observed final states with various lepton flavour contents will indicate the 
hierarchy among the A^ 33 s for different z's. 
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